home *** CD-ROM | disk | FTP | other *** search
- #LAPLACE
- #USEFLOAT=FALSE
- #SIMPLIFY=TRUE
- #ERRBND=15
- #EXPDSP=5
- #PREC=2
- #F1=f(x)=
- #F2=g(x)=
- #F3=h(x)=f(x)/g(x)
- #F4=matrix(
- #F5=vector(
- #F6=diff(f)
- #F7=diff(g)
- #F8=diff(h)
- #F9=taylor(f,0)
- #F10=taylor(f,1)
- #F11=%(run type s:startup-sequence)
- #F12=
- #F13=
- #F14=
- #F15=
- #F16=
- #F17=
- #F18=
- #F19=
- #F20=
- H=set(72.8-124,111-124,129.2-124,155-124,180.5-124,204.5-124,229.8-124,254.5-124,278-124,300.5-124,328-124)
- h_0=-256/5
- alpha=0.003674
- b=775.87
- t_z=19.38
- t(h)=(h-h_0)/(alpha*(h_0+b))
- f(t)=1+0.008*(t-t_z)/(t_z+273)
- T(h)=((h+b)*f(t(h)) - (h_0+b)*f(0))/(alpha*(h_0+b)*f(0))
- t(204)
- T(-51.2)
-
- c_w=4.182,m_w=205.3
- t_1=21.5,t_2=51
- W(t)=c_w*m_w*(t_2-t)/(t-t_1)
- dt=0.2
- W'(t)=eval(diff(W,t))
- abs(W'(49.1))*0.2
-
-